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Recently, waveguide lattices with non-uniform tunneling amplitudes have been explored due to 
their myriad tunable properties, many of which arise from the extended nature of their eigenstates 
in the absence of disorder. Here, we investigate the dynamics, localization, and parity- and time- 
reversal- (T^T) symmetry breaking in lattices that support only localized eigenstates in the disorder- 
free limit. We propose three families of tunneling profiles that lead to qualitatively different single- 
particle time evolution, and show that the effects of weak disorder contain signatures of the localized 
or extended nature of clean-lattice eigenstates. We show that in lattices with localized eigenstates, 
the signatures of T^T-symmetry breaking are acutely sensitive to the initial form of the wave packet. 
Our results suggest that waveguide lattices with localized eigenstates will exhibit a wide array of 
phenomena that are absent in traditional systems. 



I. INTRODUCTION 

Over the past decade single-mode, coupled, optical 
waveguides [l] have become a new paradigm for the real- 
ization of an ideal, one-dimensional, tight-binding lattice 
model [2j with roughly constant tunneling amplitudes 
and on-site potentials [H |5], as well as non-Hermitian 
parity- and time-reversal- (VT-) symmetric gain and loss 
potentials [6]. The electric field E{k,z) in waveguide 
/c at a distance z along the waveguide obeys equation 
of motion that is identical to that of a time-dependent 
wave function ilj{k,t); here k is the lattice-site index and 
t = z/v where v = c/ur is the speed of electromagnetic 
waves in the waveguide and nn is its refractive index. Op- 
tical waveguides have been used to simulate several phe- 
nomena such as the Anderson localization [3-5 , sponta- 
neous T^T-symmetry breaking [6 , Bloch oscillations [7 , 
quantum random walks [8l|9], and Hanbury-Brown-Twiss 
correlations due to random tunneling and on-site poten- 
tial fTO|. Most recently, one-dimensional, tight-binding 
lattice models with a position-dependent tunneling have 
been explored. The energy spectrum, density of states, 
and single-particle wave packet evolution in these tunable 
waveguides can be varied over a wide range by choos- 
ing an appropriate tunneling function [1T| Il2]; in par- 
ticular, tunneling functions that lead to commensurate 
energy levels, and the attendant wave packet reconstruc- 
tion, have been extensively explored [13, 14 . 

All of these cases have focused on lattice models that, 
in the continuum limit, describe the dynamics of a non- 
relativistic, quantum particle with a position-dependent 
mass that is determined by the tunneling function. The 
eigenstates of the corresponding Hamiltonians are, thus, 
extended [11]. Therefore, the effects of diagonal (on-site) 
and off-diagonal (tunneling) disorder in such systems, 
and their signatures in Anderson localization [15^ and 



the Hanbury-Brown-Twiss correlations [TB^, correspond 
to those in condensed matter systems. Since these sys- 
tems are typically probed through transport, tradition- 
ally, lattice models that have only localized eigenstates 
in the clean limit have not been studied; such models 
would result in a vanishing electrical conductivity even 
in the disorder- free limit. Optical waveguides, on the 
other hand, are probed via light intensity measurements 
along the waveguide, not by the motion of a wave packet 
across them. They offer an unprecedented ability to tai- 
lor the site-dependent tunneling [13, 17 ; something that 
is difficult to do in condensed matter systems and, to a 
lesser degree, in optical lattices [18", T9^. 

In this paper, we investigate the dynamics, disorder ef- 
fects and T^T-symmetry breaking in non-uniform, tight- 
binding lattice models in which most of the eigenstates 
are localized in the absence of disorder. We propose three 
classes of tunneling functions that demonstrate the wide 
range of properties typical of such models. Our three 
primary results are as follows: (i) The energy spectra 
and wave packet time-evolution have properties with no 
counterparts in systems with extended eigenstates. (ii) 
The wave packet localization due to a weak disorder en- 
codes the (extended or localized) nature of eigenstates 
even for lattices with identical clean-limit energy spectra, 
(iii) Singatures of T^T-symmetry breaking are acutely 
sensitive to the initial form of the wave packet. Our 
results show that light propagation can be significantly 
controlled and manipulated in coupled waveguides with 
localized eigenstates. 

The plan of the paper is as follows. In the following 
paragraphs, we introduce the tight binding model for a 
non-uniform N-site lattice. Section [IT] presents the prop- 
erties of a clean lattice for the three classes of tunneling 
functions that we introduce. We consider the effects of 
a weak disorder on the time-evolution in such lattices 
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in Sec. Ill In Sec. |IV[ we expand these models to their 
T^T-symmetric counterparts, and consider the effect of a 
single pair of balanced loss and gain impurities at mirror- 
symmetric positions. We conclude the paper with a brief 
discussion in Sec. [Vl 
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A lattice of coupled waveguides is described by the 
following Hamiltonian, 



N-l 



N 



where \k) represents the state with a (single) particle at 
site /c, > is tunneling amplitude between sites k 
and + 1, and Vk is the on-site potential determined by 
the local index of refraction at site k. Hamiltonian 
represents open boundary conditions, t(0) = = t{N 
Note that although t{k),Vk have the units of energy, we 
will use their scaled versions, t{k)/hc^Vk/hc^ which have 
the units of inverse- length, in the discussion of possible 
sample parameters [H [9l [17] . 

We choose tunneling amplitude profiles characterized 
by two continuous, distinct, functions. to(2/c — 1) on the 
odd sites and t£;(2/c) on the even sites are chosen such 
that tE{2k)/to{2k - 1) <C 1 for most k - 0{N). In the 
limit tE = 0^ the system decouples into N/2 waveguide- 
pairs or dimers, the energy spectrum of Eq.Q is given 
by Tto{^)^ Tto(3), . . and the eigenfunctions are local- 
ized, symmetric {S) and antisymmetric (A), dimer states 
given by \S{A),2k - 1) = {\2k - I) ± \2k))/V2. Thus, 
the particle-hole symmetric energy spectrum of the clean- 
lattice Hamiltonian is determined by the odd-tunneling 
function to- 

When tE ^ 0^ the Hamiltonian in the dimer basis be- 
comes a block-tridiagonal, symmetric matrix with the di- 
agonal block D2k-i and off-diagonal block T2k-i,2k-\-i = 



^2fe+l,2fe-l 



given by 



D2k-1 
^2/c-l,2/c+l 



to(2fe-l)diag [-1,+!]^ 

+1 +1 
-1 -1 



tE{2k) 



(2) 
(3) 



In the following section, we explore the consequences of 
such tunneling profiles for a disorder- free lattice. 



II. DYNAMICS IN A CLEAN LATTICE 

We start with three tunneling functions, 

fc2 odd, 
k even, 



tp{k) 
ti{k) 
te{k) 



k odd, 
/? ln(A;) even, 

[k{N-k)]^^^ odd, 
j3\n[k{N - k)] even. 



(4) 
(5) 
(6) 



each of which satisfies the criterion tEjto <C 1 and the 
energy-scale prefactor t' in Eqs. (|4|-([6| has been set to 
unity, t' = 1] recall that for traditional waveguide lattices 
with constant tunneling, t' /he - 10^ - 10"^ m-^ [4ll9l[T7j. 
As we will show below, the functional forms - the power- 
law tunneling function tp{k)^ logarithmic tunneling func- 
tion ti{k)^ and the tunneling function te{k) that supports 



extended states - of the tunneling amplitude broadly dic- 
tate the results that follow. 
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FIG. 1. (color online) a) Left: Typical energy spectra for lat- 
tices with tunneling functions tp{k) (blue circles), ti{k) (black 
diamonds), and te{k) (red squares), and N = 20, /3 = 1; the 
energy eigenvalues are scaled by their maximum. For tp{k) 
(blue circles) the spectrum is quadratic, for ti{k) (black dia- 
monds) it is mostly linear, and for te{k) it is linear near the 
band edges. The primary features of the spectra are deter- 
mined by the odd-tunneling function to {2k — 1) ^ tE{2k). b) 
Right: Dimensionless eigenfunctions for ground-state (bars) 
and the center-band state (stems). The top panel shows 
ground state localized at the last dimer and the center-band 
state localized at the first dimer. The center panel shows 
a broadened ground state, localized near the last dimer, and 
the center-band state localized at the first dimer. The bottom 
panel shows an extended, Gaussian, ground state centered at 
N/2 and the center-band state localized near both edges. 

The left-hand panel in Fig. [l] shows typical, particle- 
hole symmetric [21 , energy spectra for a lattice with 
A" = 20,/3 = 1 iortp{k) (blue circles), ti{k) (black di- 
amonds), and te{k) (red squares). We use the maximum 
energy £^max = — ^min to define the energy and time 
scales; it varies as £^p,max(A^) ^ t' N'^ ^ Ei^^q^^{N) ~ t' N ~ 
Ee,Ynax{N) foT N ^ 1. Notc that the energy spectrum 
for ti{k) is exactly linear when /3 = 0, and a nonzero 
j3 <^ N creates deviations from linearity at the band 
edges [20 . On the other hand, the linearity of the spec- 
trum for te{k) arises from the exactly linear spectrum of 
T{k) = [k{N - A:)] [TIl[ll[14], and the spectrum be- 
come nonlinear when P <^ 1. Note that these results are 
true for any A^ 1. Given the wide range over which the 
tunneling amplitude [9[|T7] and the number of waveguides 
in a lattice can be varied (A" - 10 - 100) [3 H [9], the 
construction of waveguide lattices with tunneling profiles 
(|4|-([6| seems feasible. 

It follows from earlier analysis that the ground state of 
Hamiltonian ([T]) is localized near the dimer with largest 
internal tunneling, max to, whereas the eigenstate near 
the center of the band is localized near the dimer with 
smallest internal tunneling, min^o; this is true regardless 
of the number of lattice sites A^. The right-hand panels 
in Fig. [l] show these two eigenfunctions. The top panel, 
with a power-law tunneling tp{k)^ shows the ground-state 
(blue bars) and the center-band state (blue stems). The 
center panel, with a log tunneling ti{k), shows a broad 
ground state (black bars) and narrow center-band state 
(black stems); this difference between the ground-state 
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sizes is because the ratio tE{N — l)/to(iV), a measure 
of perturbation away from the dimer picture, is smaher 
for tp{k) than it is for ti{k). The bottom panel shows the 
ground state (red bars) and center-band state (red stems) 
for te{k); we remind the reader that the ground state for 
the tunnehng function T{k) is a Gaussian centered at 
site N/2 [11 . Results in Fig. [l] suggest that in spite 
of the perturbation introduced by the even-site hopping 
functions t^; in Eqs. Q and (|5|, the localized nature of 
eigenstates is preserved. 
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FIG. 2. (color online) Minimum (open symbols) and maxi- 
mum (filled symbols) inverse participation ratios R(N) as a 
function of lattice size N for the three tunneling functions 
with /3 = 1; note the logarithmic scale on both axes. For 
power- law tunneling tp{k) (blue), all states are localized, as 
are those for the log-tunneling function ti{k) (black). The 
tunneling function te{k) (red squares) has both localized and 
extended states. When /3 = 0, all eigenstates - dimers - are 
localized, and as /3 increases, the minimum R values for ti{k) 
and te{k) decrease. 

To quantify this claim, we calculate the A/'-dependence 
of inverse participation ratio R{N) for all eigenstates. 
The inverse participation ratio i?^(A^) for a normalized 
eigenstate lip) = XliLi /iK) is defined as 



N 



Ri,iN) = ^\fi 



(7) 



As ^ oo, R^{N) < 1 saturates to a nonzero value 
for a localized state and vanishes, i?^(A^) ~ with 
a > 0, for an extended state [22]. Figure [2] shows 
the minimum and maximum values of R{N) for lattices 
with = 20 — 3000 and the three tunneling functions, 
Eqs. (|4|-([6|, with /3 = 1; note the logarithmic scale on 
both axes. For power-law tunneling tp{k)^ the minimum 
(blue open circles) and maximum (blue filled circles) val- 
ues of R{N) are nonzero, almost equal to each other, 
and indicate that all localized eigenstates have approx- 
imately the same size. The nonzero minimum (black 
open diamonds) and maximum (black filled diamonds) 



R{N) values for the log-tunneling function ti{k) imply 
that the ground state, although localized, is broader than 
the center-band state (see Fig.[T]). For tunneling function 
teik)^ inverse participation ratios show the existence of 
both localized (red filled squares) and extended (red open 
squares) states. These results show that most, if not all, 
states of a disorder-free lattice with these tunneling pro- 
files are localized. 

Now, we study the time-evolution of a wave packet that 
is initially confined to one (or two) waveguides. Figure [s] 
shows the amplitude A{k^t) = \{k\il)(t))\ of the time- 
evolved wave functions |'0(t)) = exp[— zi^t/fi] |?/;(0)) for 
the three tunneling functions in a disorder-free lattice 
with N = 28 sites. The horizontal axis in each panel de- 
notes time normalized by h/E^g,^ for each tunneling func- 
tion, t/{h/Emax)' The corresponding distances along the 
waveguide, for a normalized time-range t/{E^g,^/h) = 
100, are given by ^ = ct - (10 - 10^) /N'^ - 0.01 - 1 mm 
for the power-law tunneling function tp{k) and z = ct ^ 
(10 - 10^)/A^ ^ 0.3 - 30 mm for the other two tunneling 
functions. Note that waveguides with constant tunneling 
amplitude are a few mm long [4, 9 ; thus, they correspond 
to a time-range t/(£^max/^) ~ 10^ for power-law tunnel- 
ing function and a time-range oft/ (£^max/^) 100 — 1000 
for the log-tunneling function. This is consistent with the 
observation [12] that for a non-uniform waveguide array 
with a given physical length can be used to explore short- 
time or long-time behavior based on its bandwidth. 

The top panel in Fig. |3] shows A{k^ t) for tunneling 
function tp( A:), and |V^(0)) = (|7V/4) + |3A^/4))/V2. Since 
all states of this Hamiltonian are strongly and equally lo- 
calized, the partial wave packets located on two spatially 
separated dimers evolve independently. Each undergoes 
Rabi oscillations with the local frequency, and the ratio 

9 = 3^ of these frequencies is equal to the ratio of 
A-S energy gap 2to for the dimers at 3A/'/4 and N/A re- 
spectively. The (blue) side-wings in each case indicate 
the mixing with adjacent dimers due to the small, but 
nonzero, even-tunneling function tE{2k)/to{2k — 1) ~ 
1/k. The center panel corresponds to ti{k) with P = 0.5, 
and an initial state localized at the central waveguide, 
I '0(0)) = \N/2). Due to the linear spectrum, the wave 
packet is reconstructed after a time Ti = 27Th/AEi 
where AEi ^ 2E'^^max/A/^ is the (approximately) constant 
energy- level spacing (see Fig. [l]). The time-dependent 
width of the wave packet is due to the even-tunneling 
function tE/to ^ P ln(/c) /k > 1/k and it increases with /3 
for /3 < [20 . The bottom panel shows A{k,t) for the 
tunneling function te{k), which has both localized and 
extended eigenstates, and |'0(O)) = \N/2). We see that, 
although the wave packet remains confined to the dimer 
when P <^ 1, it spreads across all waveguides for a mod- 
erate value of /3 = 1.5 <C A^. In addition, as expected, 
it undergoes partial reconstructions that are reminiscent 
of (and due to) the perfect revival that occurs in lattices 
with tunneling function T{k) = [k{N - fc)]^/^ [12 [14]. 

Thus, we predict that the dynamics of a wave packet in 
the disorder-free lattices can be systematically controlled 
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FIG. 3. (Color online) Time-dependent amplitude A{k,t) for 
lattices with N = 28. Top panel: power-law tunneling func- 
tion tp{k) leads to Rabi oscillations with waveguide-index de- 
pendent frequency determined by the local gap 2to- Cen- 
ter panel: log-tunneling function ti{k) with /3 = 0.5 leads to 
periodic wave packet reconstruction due to its linear spec- 
trum. Note that increasing /3 leads to asymmetrical increase 
in the vertical spread of the wave packet, as is expected due 
to the higher even-tunneling amplitude at larger waveguide 
index. Bottom panel: tunneling function te{k) with /3 = 1.5 
has a partially linear spectrum and some extended states. 
For a fixed N, the maximum wave packet spread shown in 
the center and bottom panels increases monotonically with 
13. The bottom panel shows that for a moderate value of 
/3 = 1.5 ^ A/" = 28, the wave packet spreads over the entire 
array, even though for /3 = it is localized on the dimer. 



by a suitable value of P and the initial waveguide location. 



III. EFFECT OF A WEAK DISORDER 

For an infinite, one-dimensional system, an arbitrar- 
ily weak disorder exponentially localizes all states [15 . 
In a finite waveguide lattice, an arbitrarily weak disor- 
der localizes the wave packet to its initial waveguide mo- 
The localized fraction, characterized by the intensity at 
waveguide mo, saturates with time (or distance along the 
waveguide), but increases with the disorder strength vq. 
Therefore, the long-time, steady-state, disorder- averaged 
intensity I{k) = \{k\tlj{t))\l^ has a maximum at /c = mo 
and, for 1 <C mo < it symmetrically decays with ex- 
ponential tails away from mo [SjUl]. At short times, how- 
ever, the time- and waveguide-dependent intensity /(/c, t) 
is determined by the competition between the spread dic- 
tated by the disorder-free lattice spectrum and localiza- 
tion due to the disorder. 

In the following, we focus on the steady-state intensity 
I{k) for two lattices, with tunneling functions T{k) and 
ti{k) respectively. Both of them have (nearly) identical. 



linear, energy spectra in the disorder-free limit. However, 
all eigenstates of the former are extended and the latter 
are localized in the absence of disorder. 
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FIG. 4. (Color online) (a) Steady-state, disorder-averaged in- 
tensity in a waveguide lattice with A/" = 40, tunneling function 
T(/c), and initial wave packet \ip{0)) = \N/2) shows localiza- 
tion with exponential tails, (b) The time-dependent intensity 
I{k,t) for the same. The disorder strength is vo/Ema^ = 0.2. 
(c) Same as panel (a), but with tunneling function ti{k) as a 
function of /3. For small /3 < 1, the localized fraction is con- 
fined largely to the dimer (A^/2 — 1, N/2). With increasing /3, 
however, the weight at neighboring waveguides increases, the 
weight at the central maximum decreases, and the intensity 
profile approaches that in panel (a); note the scale- difference 
in the two panels, (d) The time-dependent intensity /(/c, t) for 
tunneling function ti{k) with (3 = 2 shows a larger fraction of 
the wave packet localized at the initial waveguide compared 
to its value in panel (b). 

Figure |4] shows results for a weak disorder, '^o/^max = 
0.2, in lattices with N = 40 sites, and tunneling func- 
tions T{k) (top row) and ti{k) (bottom row). We have 
verified that the results are independent of the number 
of realizations ~ 10^ for the disorder potential v^, 
and the type of disorder-potential distribution (Gaussian, 
uniform) as long the distribution as zero mean and vari- 
ance Vq. Panels (b) and (d) show typical time-dependent 
intensity I{k^t) for a state \ip{0)) = \N/2); horizontal 
(vertical) axis is the waveguide index (normalized time). 
Recall that in the absence of disorder, '^o = 0, the inten- 
sity time-evolution in the two panels should be similar 
to that in Fig. [3j When disorder is included, the local- 
ized fraction for tunneling function t/(/c), with /3 = 2, is 
greater than that for T{k). Panels (a) and (c) show cor- 
responding steady-state intensity I{k). Panel (a) shows 
that the wave packet is localized with exponential tails, 
as is expected for a lattice with extended eigenfunctions 
in the clean limit [3, 11\ Panel (b) shows qualitatively 
different intensity profiles as a function of /3; note the 
vertical-scale difference in panels (a) and (b). As /3 in- 
creases from 0.5 (blue circles) to 2.0 (red dashed line). 
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the intensity maximum at the center reduces, while the 
intensity at the neighboring waveguides increases. Even- 
tuahy, the intensity profile approaches that in panel (a). 

These results demonstrate the crucial role played 
by qualitatively different eigenstates of (approximately) 
isospectral Hamiltonians in the presence of disorder, and 
the change that the localization profile undergoes when 
going from one to the other. For Hamiltonians with tun- 
neling functions (|4|-(|6|, the continuum limit is a mas- 
sive particle with an internal, dimer, degree of freedom. 
Therefore, its disorder localization profile is different. 

IV. T^r- SYMMETRY BREAKING 

In this section, we will explore disorder-free lattices 
with localized eigenstates in the presence of a single pair 
of T^T-symmetric, gain and loss impurities ±^7 [6 . To 
this end, we consider a lattice with 2N waveguides and 
define the T^T-symmetric extension of the tunneling func- 
tion, 

( t(k) l<k<N, 
tvrik) = { tc k = N, (8) 

[ ^(2^- -k) TV + 1 < < 27V, 

where t{k) is either a power-law or log tunneling func- 
tion, Eqs. (|4|-([5|, and t^q-{N) = tc is arbitrary. Note 
that it is straightforward to modify this definition when 
the number of waveguides in an array is odd. The non- 
Hermitian, T^T-symmetric Hamiltonian for the system is 
then given by 

Hrr = H{trr) + n{\m){m\ - \m){m\) ^ F^^, (9) 

where Hitvr) is Hamiltonian (fTl) with the tunneling 
function 1^ or T{k) = [k{2N - fc"/^ = T{2N - /c), m is 
the position of the gain impurity, and fh = 2N + 1 — m 
denotes its mirror position where the loss impurity is lo- 
cated. Although H-pj- is not Hermitian, all of its eigen- 
values are real when 7 < 7c (m). The critical impurity 
strength 7c (^), in general, is a function of the tun- 
neling profile and the distance d = \m — fh\ between 
the loss and gain impurities. However, in the special 
case of nearest neighbor impurities, m = N ^ the critical 
impurity strength is given by 7c (TV) = tc and all en- 
ergy eigenvalues simultaneously become complex when 
7 > 7c(Ar) = tc [23 . We choose tc = t^ri^ - 1) = 
tpj-{N + 1), so that the tunneling amplitude at the lat- 
tice center is continuous. In the following, we focus on 
the evolution of time-dependent intensity across the VT- 
symmetry threshold, 7/7C = 1.00 ^0.01. Note that 
since the Hamiltonian i^-pr^ Eq.(|9|, is not Hermitian, 
the time evolution operator ex.p[—iHpj-t/h] is not uni- 
tary and the total intensity I{t) = Yl'k=i ^(^' ^) 
conserved. As in Sec. |III[ we consider two tunneling 
functions: the first, T^T-symmetric version of the log- 
tunneling function, Eq.([5|, with f3 = 0.5, has a linear 
spectrum and purely localized eigenstates; the second. 



T{k) = [k{2N - A:)]^/2 for a lattice with 2A^ sites, has a 
linear spectrum and only extended eigenstates. 
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FIG. 5. (Color online) Intensity I{k,t) for log (top row) 
and T{k) (bottom row) tunneling functions in a lattice with 
2N = 40 sites, gain and loss impurities z\zij at positions 
(20,21), and initial state |'0(O))/ = (|10) + \31))/V2 local- 
ized far away from them. The impurity strength is below the 
threshold, 7/70 = 0.99, for the left-hand column and above 
the threshold, 7/7C = 1.01, for the right-hand column. The 
top row shows that I{k,t) does not change significantly across 
the threshold, and the maximum intensity does not change 
from its Hermitian limit. The bottom row shows that the 
maximum intensity is higher than its 7 = limit since the 
wave packet comes across the impurities, and there is a small 
increase in the intensity as the threshold is passed. 

Figure [5] shows the intensity I(k^t) for a particle away 
from the impurities, with an initial state f = 

{\N/2) + |37V/2 + l))/V2 in lattice with 2N = 40 sites 
and nearest neighbor impurities at {m^fh) = (20,21). 
The left-hand (right-hand) column corresponds to im- 
purity strength 7/7C below (above) the T^T-symmetry 
breaking threshold. The top row shows that the time- 
dependent intensity I{k^t) for the T^T-symmetric log- 
tunneling function does not change appreciably as im- 
purity strength increases from 7/7C = 0.99 (left-hand 
panel) and 7/7C = 1.01 (right-hand panel); in addition, 
the maximum intensity is the same as its corresponding 
value in the Hermitian limit. These results are expected 
since the initial wave packet, localized away from the 
lattice center, does not "come across" the gain and loss 
impurities. The bottom row shows corresponding results 
for tunneling function T{k). We see that the intensity 
profile does not change significantly, although there is a 
minor enhancement for 7/7C = 1.01 (right-hand panel) 
compared to 7/7C = 0.99 (left-hand panel), and the max- 
imum intensity is higher than its 7 = value. These re- 
sults are due to the extended nature of eigenstates, which 
ensures that any wave packet "comes across" the gain and 
loss impurities. 
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FIG. 6. (Color online) Intensity /(/c,t) for log (top row) and 
T(k) (bottom row) tunneling functions in a lattice with 2N — 
40 sites, gain and loss impurities zLij at positions (20, 21), and 
initial state |'0(O))c = (| 15) + | 26)) /a/2 localized close to them. 
The impurity strength is below the threshold, 7/70 = 0.99, for 
the left-hand column and above the threshold, 7/7C = 1.01, 
for the right-hand column. The bottom row shows that that 
the maximum intensity is higher than its 7 = limit since 
the wave packet comes across the gain and loss impurities, 
and there is a small increase in the intensity as the threshold 
is passed; see the bottom row. Fig. [5] The top row shows 
that the intensity increases by an order of magnitude as the 
impurity strength 7 is varied from below the threshold (left- 
hand panel) to above the threshold (right-hand panel); see 
the top row. Fig. [5] 



How do these results change when the initial state is 
closer to the impurities? Figure |6] shows corresponding 
results for the same lattice, same impurity strengths, and 
the same time-range, but with initial state |?/^(0))c = 
(|3A^/4) + |57V/4 + l))/V2. The bottom row shows that 
the intensity profile does not change significantly from 
the left-hand panel, j/jc = 0-99 to the right-hand panel, 
7/7C = 1.01, although there is a minor enhancement in 
the right-hand panel (see the bottom row, Fig.[5|. This 
is because the weight of any extended eigenstate at the 
gain-impurity site is small, and therefore, the effect of en- 
countering the gain-site on a wave packet is small. The 
top row shows corresponding results for the log-tunneling 
function, where the initial wave packet locations are in- 
dicated by gray circles at (37V/4, 5A/'/4 + 1) = (15,26). 
The left-hand panel shows a moderate intensity enhance- 
ment below the threshold, 7/70 = 0.99. The right-hand 
panel shows a dramatic intensity enhancement above the 
threshold, 7/7C = 1.01; note the order-of-magnitude dif- 
ference in the intensity scale. This order-of-magnitude 
enhancement occurs because the dimer eigenstates, with 
energies near the band edges, have a strong weight at the 
gain impurity site. 



These results show that due to the localized nature 
of all eigenstates, the time-dependent intensity profile 
I{k^t) is acutely sensitive to parameters, such as P and 
the initial state, which control whether the wave packet 
"comes across" the loss and gain impurities. Thus, we 
predict that lattices with localized eigenstates provide a 
unique control over the violation of unitarity - how much 
and how rapidly does the intensity change from its value 
in the Hermitian limit - that has no counterpart in tra- 
ditional lattices with extended eigenstates. 



V. DISCUSSION 

In this paper, we have investigated the dynamics, dis- 
order effects, and T^T-symmetry breaking signatures in 
waveguide lattices that have, in the disorder- free limits 
a majority of localized eigenstates. We have presented 
three, novel tunneling profiles, Eqs. (|4|-([6|, that lead to 
such eigenstates. We have shown that the spatial spread 
of the wave packet and the frequency of Rabi oscillations 
can be controlled by the choice of the tunneling function 
and the initial position of the wave packet. We have also 
shown that the effect of weak disorder on such lattices 
is qualitatively different, and argued that, due to the lo- 
calized nature of disorder-free lattice eigenstates, the sig- 
natures of T^T-symmetry breaking in these lattices are 
acutely sensitive to the initial form of the wave packet. 

In this work, we have ignored the quart ic interaction 
term that, in the case of optical waveguides, arises from 
nonlinear susceptibility and in the continuum limit, gives 
rise to the nonlinear Schrodinger equation; this approxi- 
mation is justified at low intensities. Since the interplay 
between interactions and disorder has profound effects on 
the phenomenon of localization , it will be inter- 
esting to explore them in the present system. 

Here, we have only focused on T^T- symmetry breaking 
in even lattices with nearest neighbor impurities. Al- 
though Figs. |5] and |6] show the remarkable dependence 
of T^T-symmetry breaking signatures on the initial wave 
packet, T^T-symmetry breaking in even and odd lattices 
and the differences between them [25 is an open ques- 
tion. 

Our results are not dependent on the exact forms of 
the tunneling functions, Eqs. (|4|-(|6|. Indeed, for exam- 
ple, any power-law tunneling function (k^^k^) instead of 
(k^k'^) will give similar results \i fi ^ v. Similarly, the 
k — j3 ln(/c) tunneling function leads to a linear spectrum 
over a wide range /3 because ln(x) grows more slowly 
than any power of x. The robustness of these results 
implies that the small, ubiquitous variations in the on- 
site potential and tunneling amplitudes in experimental 
samples will not affect our predictions. The experimen- 
tal exploration of such coupled optical waveguides will 
deepen our understanding of lattice models with tunnel- 
ing profiles whose continuum limit is different from the 
traditional Schrodinger equation for a massive particle. 
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